Impact of boundaries on velocity profiles in bubble rafts 
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Under conditions of sufficiently slow flow, foams, colloids, granular matter, and various pastes have 
been observed to exhibit shear localization, i.e. regions of flow coexisting with regions of solid-like 
behavior. The details of such shear localization can vary depending on the system being studied. A 
number of the systems of interest are confined so as to be quasi-two dimensional, and an important 
issue in these systems is the role of the confining boundaries. For foams, three basic systems have 
been studied with very different boundary conditions: Hele-Shaw cells (bubbles confined between 
two solid plates); bubble rafts (a single layer of bubbles freely floating on a surface of water); and 
confined bubble rafts (bubbles confined between the surface of water below and a glass plate on 
top). Often, it is assumed that the impact of the boundaries is not significant in the "quasi-static 
limit", i.e. when externally imposed rates of strain are sufficiently smaller than internal kinematic 
relaxation times. In this paper, we directly test this assumption for rates of strain ranging from 
10~^ to 10~^ s~^. This corresponds to the quoted rate of strain that had been used in a number 
of previous experiments. It is found that the top plate dramatically alters both the velocity profile 
and the distribution of nonlinear rearrangements, even at these slow rates of strain. When a top 
is present, the fiow is localized to a narrow band near the wall, and without a top, there is fiow 
throughout the system. 

PACS numbers: 83.80.Iz,83.60.La,83.50.-v 



I. INTRODUCTION 



When systems are driven sufficiently far from equilib- 
rium, they often exhibit a series of transitions due to 
instabilities. This is particularly common in the flow of 
fluids, where instabilities occur at high flow rates. In 
contrast to this behavior, for sufficiently slow driving, 
complex fiuids have been observed to undergo a transi- 
tion from a purely fiowing state to a coexistence between 
a flowing and a solid-like state, i.e. shear localization 
Q, H 0,0, 0, IE 0, S li- In this context, we focus on 
complex fluids that are comprised of dense "droplets" 
(or particles) of one phase or material within a differ- 
ent continuous phase, such as foams, emulsions, granu- 
lar matter, and colloids. We are interested in the case 
where the droplets are sufficiently dense that there ex- 
ists a critical value of applied stress, the yield stress, 
below which the material does not fiow at all. In this 
situation, it has been observed that under conditions of 
non-uniform stress the material segregates into a region 
that fiows (above the yield stress) and a region that does 
not fiow (below the yield stress) [l^. However, because 
most of these materials are optically opaque, it is only 
recently that the spatial dependence of the average ve- 
locity of the "droplets" in these materials has been mea- 
sured quantitatively. For three dimensional systems, a 
key development for such studies has been the develop- 
ment of magnetic-resonance-imaging 5] techniques that 
allow for spatially resolved velocity profiles. Equally use- 
ful has been the use of quasi-two dimensional systems 
in which all the droplets can be imaged 4, 6]. Coupled 
with the experimental advances, there have been a num- 
ber of simulations that explicitly look at the possibility 
of shear localization within the context of various models 



of granular matter and foams Q • 

A striking feature of the experimental studies of shear 
localization in complex fiuids is the division of the veloc- 
ity profiles into two basic categories. The first situation 
corresponds to cases where the rate of strain is continu- 
ous across the system 0,00Q. In this case, the spatial 
dependence of the velocity is often exponential. This ap- 
pears to be the standard case for granular systems 0,0, |^ 
and bubbles confined between two plates In contrast, 
a discontinuity in the rate of strain at the transition be- 
tween the fiowing state and the jammed state is observed 
in emulsions and colloids j^, wet granular systems B , 
worm-like micelles 0, three dimensional foams [l5j |. 
and bubble rafts 0. 

In comparing the systems mentioned above, it is useful 
to note that the systems were all sheared between two 
concentric cylinders. In this geometry, there is a non- 
uniform stress across the system. This might suggest 
that the localization is due to the "simple" picture that 
part of the system is above the yield stress and part of the 
system is below the yield stress. Surprisingly, there are 
a number of ways in which the experiments suggest that 
this explanation is not sufficient. For example, some of 
the systems (especially dry granular systems 3]) clearly 
exhibit density variations that impact the fiow behav- 
ior. An understanding of these variations is necessary 
for understanding the fiow localization in these cases. In 
other studies, such as with wet granular matter, there 
are strong indications that the shear localization is the 
result of a viscosity bifurcation 7] . 

In contrast to the experiments, simulations have fo- 
cused on parallel shear. In this case, a linear velocity 
transverse to the shear is expected, and a nonlinear ve- 
locity profile is an indication of some type of shear lo- 
calization. As with the experiments, simulations exhibit 
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different behaviors depending on the details of the model. 
For example, shear localization is observed below a crit- 
ical rate of strain and under different conditions, 
above a critical rate of strain [T^ . 

For foam the situation is particularly interesting. For 
three-dimensional foam, both localized flow and flow 
throughout the system Tg', Tt] have been observed. For 
quasi- two dimensional experiments, dramatically differ- 
ent types of flow localization has been observed depend- 
ing on whether or not the bubbles were confined between 
two plates 4] or a bubble raft was used 6] . These last two 
experiments highlight the need for a systematic study of 
the impact of the confining plates when studying quasi- 
two dimensional systems. In these systems, there is al- 
ways a lower boundary supporting and confining the sys- 
tem, and depending on the experiment, there is often an 
upper boundary. Typically, the external shear is gener- 
ated by motion of the sides, with the upper and lower 
boundaries held fixed. Because the focus is understand- 
ing the behavior under conditions of small applied rates 
of strain, the systems are often described as being in a 
quasi-static limit. If true, the expectation is that the 
interaction with the confining boundaries is irrelevant. 
However, the previous experiments 0, indicate that 
the boundaries play a critical role, and suggest that one 
is not truly in a quasi-static regime, even though the be- 
havior is rate independent 18]. 

The flow behavior reported on in Refs. 0, used a 
Couette geometry, i.e. flow between concentric cylinders. 
For bubbles confined between two plates, the shear local- 
ization corresponds to an exponentially decaying velocity 
as a function of the distance from the inner cylinder 4]. 
For the case of a bubble raft (a single layer of bubbles 
floating on the surface of water [iM Ull ) ^ the velocity 
as a function of distance from the inner cylinder exhibited 
a discontinuity in the rate of strain 6] . For the case of the 
confined bubbles, simulations suggest that nonlinear re- 
arrangements of bubbles (known as Tl events) provided 
a focusing of the stress field that produced the shear lo- 
calization [l^. A Tl event corresponds to a neighbor 
switching where two neighboring bubbles separate, and 
two bubbles that were not neighbors become neighbors 
(see Fig. 1). For the bubble raft, the distribution of 
Tl events were studied and no localization was observed 

As discussed, the most striking difference between the 
two experiments is the boundary conditions on the "top" 
and "bottom" of the bubbles. The experiments in the 
confined geometry have a glass plate in contact with the 
bubbles both on the top and the bottom. For the bubble 
raft, the top surface is free, as the bubble float on a water 
surface. There is a third geometry that has commonly 
been used to study quasi-two dimensional foam: a bub- 
ble raft with a top plate in contact with the bubbles. For 
example, this has been used to study quasi- static strains 
23] and the flow around obstacles [sj. In this paper, 
we report on experimental studies aimed at determining 
the impact of the various boundary conditions. For the 



FIG. 1: (a) Schematic representation of a Tl event illus- 
trating two neighboring bubbles (A and B) switching to next 
nearest neighbor, and two next nearest neighbors (C and D) 
becoming neighbors. (b),(c),(d) A sequence of images that 
have been thresholded to bubbles are black taken from a bub- 
ble raft illustrating illustrating a Tl event. 

purposes of this comparison, we have focused on rela- 
tively monodisperse systems subjected to parallel shear. 
Monodisperse bubbles were used because these systems 
were the most reproducible between the two geometries. 
To allow for minimal variation between the systems while 
varying the boundary conditions, we focused on the two 
bubble raft systems: with and without a top. For com- 
parison with past experiments, we consider a range of 
rate of strain that was consistent with the rates of strain 
used in Ref. 0,111. 

The remainder of the paper is organized as follows. 
Section II describes the apparatus and methods for pro- 
ducing the bubble rafts in detail. Section III describes 
the method for analyzing the bubble dynamics, especially 
the identification of Tl events. Finally, Sec. IV presents 
the results and the discussion of the results. 



II. EXPERIMENTAL DETAILS 

The experimental setup contains three parts: the 
trough, the driving system and the imaging system. A 
schematic of the trough is given in Fig. [21 The trough 
consists of a rectangular Delarin dish (indicated by (A) 
in the Fig. ^ that is 300 mm x 400 mm x 75 mm. This 
serves as the main reservoir for the aqueous solution. In- 
side this dish is a Teflon frame (indicated by (B)) that 
is used to establish a symmetric boundary and can sup- 
port a glass top (not shown in Fig.I2|). The frame is held 
by four poles and controlled by four micrometers out- 
side of the trough (not shown). The size of the frame 
is 180 mm x 300 mm x 10 mm. The frame can move 
in three dimensions, with adjustments in the plane used 
to maintain symmetric lateral boundaries. The vertical 
adjustment of the frame controls the height of the glass 
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FIG. 2: (a) A schematic of the apparatus as viewed from 
the top. The details are described in the text. Highhghted in 
the figure are the driving bands (C) that are used to generate 
flow. A close up photograph of a driving band is given in (b) . 
The spacing for the bands is 4 mm. 



top relative to the bubbles. 

The bubbles are driven by two counter-rotating belts 
(indicated by (C) in Fig. [SJ using a stepper motor. As 
indicated in Fig.[2j we define the direction parallel to the 
belts to be the x-direction and the direction perpendic- 
ular to the belts as the y-direction. The stepper motor 
is a Mdrive 23 motor from Intelligent Motion System, 
Inc, model number MDMF2222, with microstepping ca- 
pability. For driving the foam, the motor is set to 51200 
microstep/rotation. The shafts, gears and belts are from 
W.M. Berg Inc. The driving bands are 210 mm long 
and spaced 57 mm apart. All the shafts are mounted at 
the bottom of the trough. The shafts are arranged such 
that the bands may be driven from outside the trough 
(through the connection indicated as (D) in Fig.[2j. This 
allows for both placement of the glass top and imaging 
the system from above the glass plate. 

The bands that are used to drive the flow act as paral- 
lel walls moving at a constant speed. They are configured 
to move in opposite directions, ensuring a location (or re- 
gion) of zero velocity in the flowing bubbles. To achieve a 
no-slip boundary conditions, belts with a groove spacing 
on the order of the average bubble size were used. The 
top of the belts are set at a height such that a single row 
of bubbles fits into the grooves on the belt. 

For imaging the system, a standard CCD camera with 
a telephoto lens is used. The lens has a focal length of 
6 mm. The focus and the aperture are manually adjusted 
to optimize image quality by minimizing distortion and 
balancing the field of view with magnification of the bub- 
bles. Images from the camera are directly digitized to the 
computer using a National Instruments frame grabber at 
a maximum frame rate of 30 frames/s. The actual frame 
rate was chosen based on the rate of strain to ensure the 
ability to track bubble motions. Selecting a frame rate 
that corresponded to a total applied strain of 0.001 be- 
tween images was found to be adequate to track bubbles 
without an excessive overload on the number of images 
required to analyze sufficiently long total strains. This 
requirement combined with the rate of strain determined 
the frame rate for any given set of images. 

The manufacture of the bubble rafts without a top is 
discussed in detail in Ref. Essentially, a solution 

of 80% of DI water, 15% of glycerin and 5% of a com- 
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FIG. 3: (a) Image of a typical set of bubbles in the case with 
no top on the system. The scale bar represents 7 mm. (b) 
Image of a typical set of bubbles for the case of a top on the 
system. The scale bar represents 7 mm. 



mercially available bubble solution ("Miracle bubbles" 
from Imperial Toy Corporation) by volume is used. Com- 
pressed nitrogen gas is flowed through the solution, with 
the flow rate and needle diameter controlling the size of 
the bubbles. By fixing the flow rate, we were able to 
generate essentially monodisperse systems. Without a 
top, the average diameter of the bubbles was 2.69 mm, 
with a standard deviation of 0.09 mm based on fitting 
the bubble size distribution to a Gaussian. The bubble 
raft is stable for about two hours without significant pop- 
ping. For producing the system with a top, the following 
procedure was used. A top plate, made from a 2 mm 
thick glass, is cleaned with a soap/water solution. Then, 
the glass is rinsed thoroughly with the same solution that 
constitutes the bubble raft in order to minimize the influ- 
ence of any transient wetting or pinning dynamics. The 
top glass is placed on the Teflon frame, completely seal- 
ing the system. When making the bubbles, the top plate 
is moved to one end of the frame, creating a small open- 
ing on the other end. Bubbles are formed at the closed 
end, driving the bubbles towards the opening. When 
the bubbles fill the entire frame, the glass top is moved 
back into position to seal the system. Again, we used a 
monodisperse system. With a top, the average diameter 
of the bubbles was 2.43 mm, with a standard deviation 
of 0.08 mm based on fitting the bubble size distribution 
to a Gaussian. 

Figure shows a typical arrangement of the bubbles 
with and without the top. Both images include the bands 
used to drive the bubbles. One can see that the region 
outside the bands is filled with bubbles as well. One 
other difference between the two systems is the nature 
of bubbles well outside the bands. For the case of a top, 
because the system is effectively sealed, the bubbles fill 
the entire region within the supporting frame. For the 
case without a top, the flow outside the bands does show 
some unavoidable multilayer formation. This occurs in 
the corners of the Teflon frames. The loss of bubbles to 
these multilayers results in the formation of voids on the 
inner perimeter of the Teflon barriers. The density of 
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bubbles between the bands, in the region of interest, is 
not noticeably affected by this. While the multilayers and 
void formation may have consequences for the pressure 
or stress fields globally, we find our velocity profiles and 
Tl rates/densities do not depend on the occurrence or 
growth of the multilayers or voids. We have checked for 
variation along the x-direction in many of the system 
properties due to the influence of the bubbles outside the 
flowing region. We observe a small entrance effect that 
decays rapidly. Therefore, we focus on the central area 
of the driven region. 



III. ANALYSIS METHODS 

The primary dynamical features of the system we ex- 
tract are the velocities of the individual bubbles and local 
topological rearrangements. The main topological event 
of concern for this paper are the Tl events, where neigh- 
bor rearrangements occur (see Fig. 1). 

The raw data from the experiments consists of an im- 
age series capturing the time evolution of the bubbles at 
different rates of strain. The analysis of these images may 
be classified in two sections (1) A reduction of each image 
to a set of bubble centers, edges and vertices and (2) The 
evolution of these reduced measures between successive 
images to extract velocity profiles and Tl events. 

The images are initially cropped to a desired region 
of interest and Fourier filtered to eliminate noise associ- 
ated with the CCD camera and optical non-uniformities. 
The grayscale images are then reduced to binary images 
by thresholding them at an appropriate value to demar- 
cate the interior regions of bubbles/cells from the bubble 
edges. The positions of the centers of each bubble are 
computed as the centers of mass of the interior regions of 
the cells in such a binary representation. This procedure 
reliably identifies over 99% of the bubbles in each image. 

The center positions in consecutive images that show 
the least displacement are identified as being associated 
with the same bubbles. To reliably make such identifi- 
cation requires the displacement of the bubbles between 
successive images be less than their radii. This was one 
criteria used in selecting the frame rate. The velocity 
of the bubbles is computed using the displacement of 
the bubbles between two images and the time taken for 
the displacement and averaging over many bubbles and 
frames. 

For the purposes of this paper, the velocity profiles 
represent an average over a total applied strain of 5 and 
a spatial average in the x-direction. The velocity profile 
is essentially independent of the x-position in most of the 
central region of the trough. There is a small entrance 
length at each end in which the velocity profile varies. 
Therefore, to be conservative, only the central 1/3 of the 
trough (in the x-direction) is used for computing aver- 
age velocities. To confirm whether or not slip exists at 
the driving bands, we computed velocities for the entire 
width of the trough (in the y-direction). The y-direction 




FIG. 4: (color online) Illustration of identifying Tl events 
using the Voronoi construction. Image (a) and image (b) are 
the Voronoi reconstruction for two sequential images. Bubbles 
involved in Tl events are colored. 



is divided into evenly spaced bins, and all bubbles in a 
given y-bin, independent of their x position, are used to 
compute the average velocity at that point. 

In our experiments, the the bubbles form densely 
packed two dimensional structures. We build a space fill- 
ing tessellation from the positions of the cell centers us- 
ing a Voronoi construction. The edges and vertices thus 
extracted are seen to accurately reproduce the network 
formed by the edges of individual cells in the bulk of the 
system. At the boundaries of the network, the Voronoi 
construction is not representative of the cell edges and 
therefore in any further analysis we expunge cells that 
have vertices at the boundaries. 

Knowing the vertices shared between cells makes it 
possible to identify cells that are neighbors of each other. 
In order to identify Tl events occurring in the system, 
we identify cells for which the next-nearest neighbors be- 
come nearest neighbors. This scheme identifies two of 
the cells that participate in a Tl-event. The other two 
cells correspond to those in which nearest neighbors be- 
come next-nearest neighbors. While this methodology 
of identifies pairs of cells participating in Tl events, a 
number of such pairs often occur in proximity forming 
clusters that may be associated with slip zones. The size 
of these clusters seen depends on the framerate of image 
capture. However, assuming one had a sufficiently fast 
camera, all individual Tl events might be observed. The 
positions of the Tl events may be computed as the cen- 
ter of mass of the cells in each cluster. An example of 
the Voronoi reconstruction and detection of Tl events is 
given in Fig. 01 The two images illustrate the system be- 
fore and after Tl events occur. Bubbles involved in the 
Tl events are shaded for easy identification. 

The Tl events correspond to regions where slips be- 
tween cells occur resulting in neighbor switching. These 
events are the primary mechanisms through which flows 
in foam systems are known to occur. The Tl events re- 
flect a variation in the connectivity between neighboring 
cells from as a metric independent measure, while the 
velocity profiles of the bubbles are based on a eucledian 
metric. The relationship between the externally imposed 
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FIG. 5: (color online) The scaled velocity profile {U = 
(bubble velocity) /(belt velocity)) as a function of scaled po- 
sition {y/ < d >) across the trough for three different rates of 
strain for a system without a top, where < d > is the average 
bubble diameter. The scale of the position axis is set so that 
it extends from the edge of one belt to the other belt. 



shear inducing local Tl events and velocity profiles are 
explored in the next section. 



IV. RESULTS 

The main result of the paper is a comparison of 
Fig. [51 and Fig. [SI In both figures, three different 
scaled velocity profiles are plotted as a function of the 
displacement from the center of the trough in the y- 
direction (scaled by the average bubble diameter). (The 
velocity is scaled by the driving belt velocity, U = 
(bubble velocity) /(belt velocity).) A few common fea- 
tures of the velocity profiles are worth highlighting. If 
there is no slip at the boundary, the scaled velocity should 
be one by definition. Second, the velocities scale for 
both boundary conditions and the three rates of strain 
reported on here. This indicates that we are in a rate 
independent regime Finally, because the bands are 
moving in opposite directions, the velocity goes through 
zero, and it is expected to be zero in the center of the 
trough. Both profiles are consistent with this expecta- 
tion. 

The most striking feature is the extreme localization 
of the flow when there is a top (Fig. and the corre- 
sponding essentially linear profile without a top (Fig.jS}. 
This provides strong evidence for the importance of ac- 
counting for the confining boundary conditions, even in a 
case where one expects the rate of strain to be sufficiently 
slow. For comparison with earlier work, the data with a 
top is plotted semi-log in Fig. [7| One can see that the 
behavior near each boundary is consistent with an expo- 



FIG. 6: (color online) The scaled velocity profile {U = 
(bubble velocity) /(belt velocity)) as a function of scaled po- 
sition {y/ < d >) across the trough for three different rates 
of strain for a system with a top, where < d > is the average 
bubble diameter. The scale of the position axis is set so that 
it extends from the edge of one belt to the other belt. 



nential decay over a few bubble diameters. The devia- 
tions from the exponential behavior in the center may be 
in part due to the experimental resolution of our velocity 
measurements. Also, it should be noted that the profile 
in the case with no top (Fig.[5J is not perfectly linear, as 
would be expected for a "simple" fluid. One candidate for 
the deviations from linearity is the monodispersity of the 
bubble raft. This is certainly an interesting question, and 
will be the subject of future more detailed work. How- 
ever, for the purposes of establishing the impact of the 
boundaries, the difference between the profiles in Fig. [51 
and Fig. El are more important than the variations from 
linear velocity in the case of not having a top. 

The other feature of the flow that is apparent is the 
behavior at the driving band. In the case of no top, we 
achieved a no-slip boundary condition by containing the 
bubbles in the spaces in the bands. We tested this by 
varying the position of the band relative to the bubbles. 
If the height of the band was such that the bubbles sat 
at the edge of a band but not in one of the gaps, we 
observed complete slip at the boundary. In this case, no 
flow was observed anywhere in the system. For the case 
with the top, we observed some slip at the boundary. 
However, the "slip" was not complete in the sense that 
the bands were still able to drive flow, just with a reduced 
average speed relative to the speed of the bands. The 
introduction of slip in the case of the top is most likely 
the result of the drag from the top acting on the bubbles. 
An interesting feature of the slip is that the degree of slip 
was independent of the rate of strain. This suggests that 
both the force between the bubbles and the driving band 
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FIG. 7: The natural log of the scaled velocity profile (ln([/)) 
as a function of scaled position {y/ < d >) across the trough 
for two rates of strain in the system with a top (solid squares 
are a rate of strain of 0.014 s~^ and open triangles are 



0.0014 s" 



The lines represent linear fits to the 0.014 s 



data, indicating an exponential decay of the velocity in that 
regime. 
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FIG. 8: Positions of Tl events in the central portion of the 
system (a fraction of the x direction) without a confining top. 
Both the X and y position is scaled by the average bubble 
diameter < d >. 



and the drag of the plate on the bubbles are independent 
of rate of strain. One could test the impact of the plate in 
the future by selecting driving bands with varying degrees 
of interaction between the band and the bubbles. For a 
sufficiently strong interaction, one would expect no slip, 
despite the drag due to the plate. 

To further explore the impact of the confining top 
boundary. Fig. [SI and El compare the spatial distribution 



FIG. 9: Positions of Tl events in the central portion of the 
system (a fraction of the x direction) without a confining top. 
Both the X and y position is scaled by the average bubble 
diameter < d >. 



of Tl events. For these plots, only a total applied strain 
of 0.5 is used. (The smaller interval of strain is used to 
avoid overcrowding the plot.) To ensure that the steady 
state statistics are being viewed, the last 0.5 of strain 
out of a total strain of 5 is selected. Because we are 
interested in the differences of the boundary conditions 
at the slowest possible rate of strain, only the case for 
7 = 1.4 X 10~^ s~^ is shown. Each circle represents the 
spatial location of a Tl event. 

The most dramatic feature is the absence of Tl events 
from the center region when a top is placed on the sys- 
tem. This is highlighted in Fig. 11771 where the probability 
of a Tl event occurring at a particular y/ < d> position 
in the range —22<x/ < d > < 22. The histogram illus- 
trates the dramatic difference between the number of Tl 
events in the central portion of the system for the two 
cases. This indicates the strong connection between the 
occurrence of Tl events and the existence of non-zero ve- 
locity. Also of interest, is the slight increase of Tl events 
near the boundaries that results in a peak in the proba- 
bility for both cases. The peak is close to the boundary, 
and more pronounced for the case with the top. The 
nature of the peaks reinforces the relative boundary slip 
for the two cases, as the probability of Tl events right 
at the boundary drops to zero for the case without a 
top, confirming a lack of slip. In contrast, with the top, 
some percentage of Tl events occur even very close to 
the boundary, as must happen if slip occurs. 

A noteworthy feature of the distribution of Tl events 
is the intermittent occurrence of coherent events along 
lines throughout the system. These structures appear 
as lines in Fig. [SI and correspond to relaxation through 
a number of neighbor rearrangements resulting in large 
scale slip zones within the bubble raft. (Various views 
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FIG. 10: (Color online) Histogram of the probability of a 
Tl event (P(T1 event)) as a function of ^/ < d > computed 
using in the range —22 < x/ < d > < 22 for both the case 
with a top (solid red bars) and without a top (hashed open 
bars). The histogram shows the absence of Tl events from 
the center of the system with a top present, and the relative 
flatness of the distribution in the center for flow without a 
top. 



of a three-dimensional space-time plot of the occurrence 
of Tl events is available on EPAPS jl^ that illustrates 
the temporal correlations between events.) Preliminary 
results with a high-speed camera suggests that the chains 
of Tl events align with the crystallographic axes of the 
bubble lattice, but more systematic work is necessary to 
confirm this correlation. Also, it is interesting to spec- 
ulate on the correlation between these apparent "slip 
planes" in which the majority of Tl events occur and 
the observed systematic variation in the velocity profile 
from linear (see Fig.ISJ. Clearly, future work is needed to 
both clarify the impact of monodispersity on the degree 
of linearity of the velocity profile and to explore more 
quantitatively the connection between Tl events and ve- 
locity profiles. 

In summary, our results demonstrate that for otherwise 
identical systems under conditions of small applied rates 
of strain, the existence of a confining solid plate can have 
a dramatic impact on the flow behavior of the system by 
suppressing flow in most of the system. The suppression 
of flow is paralleled by a suppression of Tl events, con- 
firming the expected strong connection between the two 



processes. This strongly suggests that even when there 
is evidence of rate independence (for example, the scal- 
ing of velocity profiles as seen in this paper), one has to 
be very careful in interpreting the measurements from a 
system with a top in place. 

In trying to understand the role of the top plate, the 
results point to three primary sources of damping for 
the bubble rafts considered. These correspond to viscous 
drag between (a) bubbles, (b) bubbles and the water sub- 
phase and (c) bubbles and the glass top plate (Fig. or 
air (Fig. 0). The bubble-bubble interaction can be con- 
sidered to be the "intrinsic" dissipation that provides the 
effective viscosity for the system under fiow. Previous ex- 
periments indicate that the bubble-bubble interactions 
are significantly stronger than any bubble-subphase in- 
teraction The velocity profiles indicate the relative 
strength between (a) and (c). With a top, the exponen- 
tial decay of velocity at the boundaries (Fig[7I) suggests 
that the viscous damping between the bubbles and the 
top glass plate dominates relative to the bubble-bubble 
interactions, resulting in an exponential profile. Under 
the current conditions, without a top, the bubble-bubble 
interactions produce a velocity profile that is consistent 
with a model that treats the foam as a viscous fluid, even 
if the viscosity is non-Newtonian. What remains to be 
seen is under what, if any, conditions the velocity profile 
in the linear case exhibits significant departures from lin- 
ear. For example, at sufficiently slow rates of strain or for 
larger system sizes, new behavior may be observed. Fi- 
nally, it should be noted that recent work on a system of 
only two bubbles connects these differences to whether or 
not a system is modelled as a dry foam (confining plates) 
or a wet foam (bubble raft without a top) 27], which can 
also be connected with the different mechanisms of dis- 
sipation. This is an interesting connection that requires 
further exploration in the large system. 
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